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Abstract 

The problem of finding low rank m x m matrices in a real affine subspace of 
dimension n has many applications in information and systems theory, where low 
rank is synonymous of structure and parcimony. We design a symbolic computation 
algorithm to solve this problem efficiently, exactly and rigorously: the input are the 
rational coefficients of the matrices spanning the affine subspace as well as the 
expected maximum rank, and the output is a rational parametrization encoding a 
finite set of points that intersects each connected component of the low rank real 
algebraic set. The complexity of our algorithm is studied thoroughly. It is essentially 
polynomial in where r is the expected maximum rank; it improves on 

the state-of-the-art in the field. Moreover, computer experiments show the practical 
efficiency of our approach. 
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1 Introduction 

1.1 Problem statement 

Let Q, M and C be respectively the fields of rational, real and complex numbers. Let m, 
n, r be positive integers with 0 < r < m — 1 and let do,..., be m x m matrices with 
entries in Q. Let x = {xi,, Xn) be a set of n variables. We consider the affine map (or 
linear matrix) A[x) defined by 


(Xi, . . . , Xn) I—t A[x) — do -|- XiAi -!-■■■+ XnAn. 
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Given A{x) as above, we consider the following algebraic set: 

Vr = {x e \ rank A{x) < r}. 

The goal of this paper is to design an efficient exact algorithm for solving the problem 
of computing at least one point in each connected component of the real algebraic set 
Vr n M”. 

By exact algorithm, we mean that we rely on multi-precision arithmetic; our output 
which encodes points whose coordinates are algebraic numbers is exact since it provides 
a rational parametrization with rational coefficients of those points. However, we allow 
probabilistic algorithms using on some random (i.e. generic) changes of variables. 

By efficient algorithm, we mean that the theoretical worst-case complexity should be well 
understood and controlled, improving the state-of-the-art in the field. In particular, our 
algorithm should perform significantly better than general-purpose symbolic computation 
available algorithms, and this should be evidenced by explicit, reproducible experiments. 

1.2 Motivations 

The problem of finding low rank elements in a given affine subspace has many applica¬ 
tions in systems, signal and information engineering, where low rank elements typically 
correspond to sparsity and structure requirements. For example, in the context of semidef- 
inite programming (SDP) hierarchies for polynomial optimization [12], low rank moment 
matrices provide guarantees of global optimality of a convex relaxation of a nonconvex 
optimization problem. 

Moreover, solving efficiently the real root finding problem for linear matrices is a first step 
towards obtaining complexity bounds for semidefinite programming. Indeed, let A{x) be 
symmetric, and let S' = {x G M" | A{x) ^ 0} be the associated spectrahedron. One can 
prove that if r is the minimum rank attained by A{x) over (the boundary of) S, then at 
least one of the connected components of 0 MA is contained in S. This implies that 
one can decide the emptiness of S (and, hence, the feasibility of a semidefinite program) 
by deciding the emptiness of Vr 0 M”. 

Similarly, the geometry of low rank structured matrices (e.g. Hurwitz matrices, Hankel 
matrices, Toeplitz matrices, resultant matrices) is pervasive in algebraic approaches to 
information engineering (including systems control, signal processing, computer vision 
and computational geometry), see e.g. [lU], [3H] or [12] and the references therein. The 
specific geometry of low rank manifolds can be exploited to design efficient nonlinear local 
optimization algorithms [1]. Sparsity-promoting optimization methods are now commonly 
used in floating-point computational environments, and compressed sensing algorithms 
based on large-scale convex optimization methods are listed amongst the success stories 
of applied mathematics in engineering, see e.g. [13]. Finally, linear matrices and their loci 
of rank defects are the object of the so-called low rank approximation problem, see e.g. 

H- 

In our paper, we are not after large-scale problem instances solved approximately with 
floating point arithmetic. In contrast, our focus is on symbolic computation and rigorous 
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algorithms. This means that we are not concerned with nnmerical scaling and condition¬ 
ing issues: all our computations are carried out with exact arithmetic on integers and 
rational numbers, and we provide mathematical guarantees of exactness of the output of 
our algorithm, under the assumption that the input is also exactly provided in rational 
arithmetic and satishes some genericity assumptions that are specihed below. Obviously, 
these guarantees come with a price, and our algorithm complexity is exponential in the 
number of variables or problem size, and hence limited to small dimensions. But this 
is not specihc to our algorithm, this limitation is shared with all symbolic computation 
methods: our algorithm should be applied to small-size problems for which it is absolutely 
crucial to hnd exact solutions. 

However, the main difference with the state-of-the-art is that the complexity achieved by 
our algorithm is essentially quadratic in a multilinear Bezout bound on the maximum 
number of complex solutions encoded by the output. This bound is itself dominated by 
• Hence, for particular sub-classes of the problem, for example when the size 
of the matrix is hxed, the multilinear bounds (and hence the complexity) are polynomial 
in the number of variables. 


1.3 Example 


Consider the Cayley determinantal cubic surface 


P 2 = I rank A{x) < 2}. 


with 


/I a:i Xa \ 

A{x) = xi 1 xa 

\ X2 X3 1 / 


whose real part delimits a convex open connected component whose closure is the well- 
known spectrahedron arising in the SDP relaxation of the MAXCUT combinatorial op¬ 
timization problem, see e.g. (SUl Example 2] and Figure The Cayley surface has 
singularities captured by the set 


Vi = [x E \ rank A{x) < l}. 
whose real part consists of four points 


T>i n = 





1.4 State of the art 

Computing real solutions of systems of polynomial equations, and deciding the emptiness 
of real algebraic sets, is a central question in computational geometry. Since one typically 
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Figure 1: Cayley determinantal cubic surface with its four singular points on the boundary 
of its spectrahedron. 


deals with positive dimensional solution sets, one possible approach is to design algorithms 
computing a hnite set intersecting each connected component. While the complexity of the 
hrst algorithm solving this problem (Tarski, [57]) was not elementary recursive, Collins 
designed in [T6| the Cylindrical Algebraic Decomposition algorithm, whose complexity 
is doubly exponential in the number of variables. Since Thom-Milnor bound for the 
maximum number of connected components of a real algebraic set (see uni Theorem 
7.23]) is singly exponential in the number of variables, tremendous efforts have been 
made to obtain optimal complexity bounds. 


Grigoriev and Vorobjov introduced in |3l] the hrst algorithm whose complexity is singly- 
exponential in the number of variables n. It is based on the critical point method. The 
algorithms in [5]|ll|7l[6l[8l|53l|5l| also rely on the computation of critical points. On inputs 
of degree < d, they lead to almost optimal complexities which are essentially cubic in 
for the general smooth case, quartic in for the general singular one. These techniques 
have also been used in the context of polynomial optimization [331135] . 


In the context of algebraic varieties dehned by the vanishing of the determinant of a 
linear matrix, which corresponds to the problem stated in Section ]_T in the particular 
case r = m — 1, we designed in 1371 an efficient algorithm based on the study of incidence 
varieties. In this previous work, the proof of correctness of the algorithm is based on 
regularity properties of the incidence varieties, and our theoretical analysis of worst-case 
complexity also strongly exploits the determinantal structure. 


In the current paper, we build on our previous work [37[ to allow the rank of the linear 
matrix to drop below the value m — 1, which means that the algebraic variety is dehned 
not only by the vanishing of the determinant, but also that of lower order minors of the 
matrix. Consequently, the proof of correctness and the theoretical analysis of complexity 
of the algorithm should be adapted to this case. The basic ingredients of our algorithm 
remain essentially the same as those of EH: we still consider critical points of some 
well chosen projections restricted to incidence varieties. However, the arguments used 
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in 1371 cannot be applied mutatis mutandis to our context since we allow higher rank 
defects. Hence, our analysis of regularity properties of the indicence varietes as well as the 
bounds on the complexity are significantly revisited. In particular, for the computation of 
rational parametrizations we rely on homotopy-based techniques for obtaining geometric 
resolutions of finite sets, see m- 

1.5 Paper outline and main results 

The algorithm described in this paper, with input a linear matrix A{x) = Aq+xiAi + ...+ 
XnAn, with Ai G i = 0,1,... ,n, and an integer r < m — 1, computes a rational 

parametrization of a finite set intersecting each connected components of Vr fl MA. The 
design of the algorithm is intented to take advantage of the special structure of the input 
problem and hence to behave better than algorithms based on the critical point method 
that solve more general problems. 

To do that, we investigate properties of an incidence variety and show that our problem 
can be reduced to compute hnitely many critical points of the restriction of a well-chosen 
projection to this incidence variety. The system that defines these critical points has a 
special sparisty structure. Using routines in [30], one can compute a rational parametriza¬ 
tion of these critical points by exploiting this sparsity structure. We establish a bound 
6 on the degree of the parametrization, and, using |1D], we show that the complexity is 
essentially quadratic on 6. This bound is dominated by 

Moreover, we provide computer experiments that show that our strategy allows to tackle 
problems that are unreachable by implementations of other generic algorithms based on 
the critical point method. 

The paper is structured as follows. In Section we set up the general notation used 
throughout the paper and we recall the key notion of incidence variety. We also state 
formally the genericity properties under which our algorithm is guaranteed to provide 
a correct output. Finally, we decribe the input/output data representation of our algo¬ 
rithm, solely based on integer and rational arithmetic. In Section we provide a formal 
description of our algorithm, we state its correctness, and we carry out a precise com¬ 
plexity analysis. The proof of correctness relies on the following technical ingredients: 
regularity of an indidence variety, see Section dimension of a Lagrangian variety, see 
Section and closure properties, see Section The paper ends up with some computer 
experiments on an implementation of our algorithm, reported in Section 


2 Definitions and notation 

We denote by Q" (resp. C") the set of vectors of length n with entries in Q (resp. C). 

A subset V C C” is an affine algebraic variety (equivalently affine algebraic set) if it 
is the common zero locus of a system of polynomials / = (/i,...,/g) G Q[x]'^, with 
X = (xi, ... ,Xn)- We also write V = /“^(O) = Z{f). Algebraic varieties in C" define the 
closed sets of the so-called Zariski topology. Zariski open subsets of C"" are sets whose 
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complement are Zariski closed; they are either empty or dense in C"". 

The set of all polynomials vanishing on an algebraic set V is an ideal and it is denoted 
by /(V) C Q[a;]. This ideal is radical (i.e. G /(V) for some integer k implies that 
g G /(V)) and it is generated by a hnite set of polynomials, say / = (/i,..., fp). We also 
write J(V) = (/i,..., /p) = (/) when a set of generators is known. We say that the length 
of the polynomial system / = (/i,..., /p) is p. 

Let GL(n, C) (resp. GL(n, Q)) be the set of non-singular n x n matrices with entries in 
C (resp. Q). The identity matrix is denoted by J„. Given a matrix M G GL(n, Q) and 
a polynomial system a; G C” h- )■ /(x) G we denote by / o M the polynomial system 
X E ^ f{Mx) G C^. If V = Z{f), the image set Z{f o M) = {x G C” : f{Mx) = 
0} = {M~^x G C” : f{x) = 0} is denoted by M~^V. Given q < n and M G we 

denote by minors(g, M) the set of determinants of g x g submatrices of M. 

For / C Q[a;]'^, we denote by Df the Jacobian matrix of /, that is the q x n matrix 
Df = When / generates a radical ideal, the co-dimension c of Z(f) is the 

maximum rank of Df evaluated at points in Z{f). Its dimension is n—c. The algebraic set 
V = Z{f) is said irreducible, if it is not the union of two algebraic sets strictly contained in 
Z{f). If V is not irreducible, it is decomposable as the hnite union of irreducible algebraic 
sets, called the irreducible components. If all the irreducible components have the same 
dimension, V is equidimensional. The dimension of V coincides with the maximum of the 
dimensions of its components. 

Let /: C” —)■ generate a radical ideal, and let V = Z{f) be equidimensional of di¬ 
mension d. A point a; G V such that the rank of Df is equal to n — d is a regular point, 
otherwise is a singular point. We denote by regV and singV respectively the set of regular 
and singular points of V. 

Let /: C” —)■ and let V = Z{f) be equidimensional of dimension d. Let g: C” —)■ 

and assume that / generates a radical ideal. The set of critical points of the restriction 
of g to V is dehned as the set of points at which / and the minors of size n — d + m of 
the extended Jacobian matrix D{f,g) and at which D{f) has rank n — d. It is denoted 
by crit(g,V). Let tti : C” —)■ C be the projection tii{x) = xi and let Dif be the matrix 
obtained by deleting the hrst column of Df. Then crit(7ri, V) is equivalently dehned by 
the zero set of / and by the maximal minors of Dif. 


2.1 Incidence variety 


Let A = (Aq, Ai, ..., An) he m x m matrices with entries in Q, and A{x) = Aq -f- XiAi + 
■ ■ ■ + XnAn the associated linear matrix. \f x E = {x E \ rank A(a;) < r}, the kernel 
of A{x) has dimension >m — r. We introduce m{m — r) variables y = (gi,i,..., ym,m-r), 
stored in a m x {m — r) linear matrix 


y{y) = 


Vip ' ' ' yi,m—r 


Vm,! ' ' ' ym,m—r J 
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and, for U G and S G r)x{m. r) define the incidence variety associated 

to (A, U, S) as 

VM,U,S) := {{x,y) ■. A{x)Y{y) = t),UY{y) - S = t)]. (1) 

Remark 1 It is easy to check that the matrix Y[y) has full rank m — r if and only if 
there exists U G of full rank and S G GL(m — r, Q) such that UY (y) — S = 0. 

For A G U = {Uij)i<i<m-r,l<j<m € Q("^-’')x™ and S = e 

Q(m-r)x(m-r)^ define 

U Q'j ■ (^n+m(m—r) _ y f^(2m—r)(m—r) 

(x,y) ^ (A(x)Y(y), [/Y(y)-S) 

Remark that Vr{A, f/, S) = Z{f{A, U, S)) and that the projection of Vr{A, U, S) over the 
x—space is contained in the determinantal variety Vr, by definition. We will prove that 
np to genericity assnmptions, the algebraic variety Vr{A,U, S) is eqnidimensional and 
smooth. When parameters are clear from the context, we will denote f{A, U, S) by /. 


2.2 Data representation 

2.2.1 Input 

The inpnt of onr algorithm is a m x m linear matrix A{x) = Aq + xiAi + ■ ■ • + XnAn 
encoded by the vector of defining matrices (Aq, ^i,..., An), with coefficients in Q, and an 
integer r snch that 0 < r < m — 1. The vector (Aq, Ai,..., An) is nnderstood as a point 
in . With abnse of notation, we refer to this inpnt with the short-hand notation 

A = (^0; Ai,..., An). 


2.2.2 Output 


The output of the algorithm encodes a finite set intersecting each connected component of 
Indeed, the initial problem is reduced to isolating the real solutions of an algebraic 
set ^ C C"" of dimension at most 0. To this end, we compute a rational parametrization 
of Z that is given by a polynomial system q = {qo,qi,... ,qn, qn+i) £ such that 

go, Qn+i are coprime (i.e. their greatest common divisor is constant), Z is in one-to-one 
correspondence with the roots of qn+i, and 





’go(fj) °}- 


This allows to reduce real root counting isolation to a univariate problem. Note that 
when qn+i is square-free, the cardinality of Z is the degree of qn+i, we denote it by deg q. 
When Z is empty, we have g = 1. 

One could remark that, since go and qn+i are coprime, one can obtain a rational parametriza¬ 
tion with go = 1. As observed in [2], choosing go equal to the derivative of g„+i leads, most 
of the time, to rational parametrizations with coefficients of smaller size in gi,..., g„. 
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Given / C Q[a;] defining a finite set ^ C C"’, there exist many algorithms for computing 
a rational parametrization of Z. 

2.3 Genericity properties 

Our algorithm works under some assumptions on the input A and r. We will prove below 
in Section that these assumptions are generic. 

2.3.1 Property Gi 

We say that A satishes Gi if singH^ = In Proposition below, we prove that this 

property is generic in the space of input matrices. 

2.3.2 Property G 2 

We say that a polynomial system / = (/i,..., /p) G Q[xi,..., Xn]^ satishes G 2 if 

• (/) is a radical ideal, and 

• Z{f) C C” is either empty or smooth and equidimensional. 

2.3.3 Property G 

Finally, we say that (A, U, S) satishes G if: 

• A satishes Gi, and 

• Vp is empty or has co-dimension [m — p)^, for 0 < p < r, and 

• f{A,U,S) satishes G 2 . 

3 Algorithm: description, correctness, complexity 

In this section, we describe the algorithm LowRank, prove its correctness and estimate its 
complexity. 


3.1 Formal description 


The input of LowRank is a couple {A,r), where A is a tuple of n -|- 1 square matrices 
Aq, Ai,An of size m, with entries in Q, and r < m — 1 is an integer. The algorithm is 
probabilistic and, upon success, its output is a rational parametrization encoding a hnite 
set of points intersecting each connected component of the real algebraic set {x G M"" | 


rankA(a;) < r}, see section 2.2 for a comprehensive description of the input and output 
data. 



3.1.1 Notation 


Change of variables. Let M E GL(n,C). As already explained in Section]^ we denote 
by A o M the affine map x h- A{M x) obtained from A by applying a change of variables 
indnced by the matrix M. In particnlar A = A o For M E GL(?7,,C), and for 
all A E g g ^ GL(m — r, C), we consequently denote by 

f{AoM, U, S) the polynomial system /(A, U, S) applied to (M x, y), and by Vr(Ao, U, S) = 
Zif{AoM,U,S)). 


Fibers. Recall that given A E M G GL(?7,, C),U E and S E GL(m — 

r, C), the polynomial system /(A o M, [/, S) and its zero locus Vr{A o M, U, S) have been 


dehned in Section 2A Given t E C, dehne 


ff ' _ y 

{x,y) ^ 


^(2m—r)(m—r)+l 

{f{AoM,U,S), x^-t) 


and denote by Vr,t(A o M, f/, S) = Z{ft) C the section of Vr with the linear 

space dehned by xi — f = 0. When parameters are clear from the context, we use 
the shorter notation Vr,t- We remark that, e.g. for M = ft = f{A,U,S) with 

A = (Aq + tAi) + X 2 A 2 + ... + XnA^. 


Lagrange systems. Given v E dehne 

i{A O M, U, S, V) : Cn+(3m-r)(m-r) ^n+(3m-r)(m-r) 

(x,y,z) H-t (f(AoM,[/,S),z'Dif,v'z-l) 

where variables z = (zi,..., Z( 2 m-r)(m-r)) stand for Lagrange multipliers, and let Z(A o 
M, U, S, v) = Z(£(A O M, U, S, v)) C C^+m{m-r)+i2m-r}{m-r) _ 


3.1.2 Subroutines 

The algorithm LowRank uses diherent subroutines, described as follows. 

• IsReg: inputs parameters A, U, S and outputs true if (A, U, S) satishes G and false 
otherwise; 

• RatPar: inputs a zero-dimensional polynomial system / and returns a rational 
parametrization of its solutions; 

• Project: inputs a rational parametrization of a hnite set Z C and a subset of 
the variables Xi,... ,X]\f, and outputs a rational parametrization of the projection 
of Z on the space generated by this subset; 

• Lift: inputs a rational parametrization of a hnite set Z C and a number t E C, 
and outputs a rational parametrization of {(f, x) : x E Z}; 
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• Image; inputs a rational parametrization of a finite set Z C C'^ and a matrix 
M G GL(iV, C); outputs a rational parametrization of M~^Z-, 

• Union; inputs rational parametrizations encoding finite sets Zi,Z 2 and outputs a 
rational parametrization of Zi U Z 2 . 


3.1.3 The algorithm 

With input {A,r), the algorithm draws randomly matrices U and S and checks whether 
A satishes G. If this is the case, it calls a recursive subroutine called LowRankRec with 
the same input {A,r). Otherwise it returns an error message. 


LowRank(A, r) 

Input: A n—variate linear matrix A{x) of size m, encoded by the m^{n+ 
1) rational entries of Aq, Ai,An, and an integer 0 < r < m — 1; 
Output: Either an error message, or the output of LowRankRec. 

Procedure: 

1. Choose randomly U G ^ Q(„j_^)x(m-r-). 

2. If lsReg(A, U, S) = false then return an error message saying that 
the input data A are not generic; 

3. else return LowRankRec(A, r). 


The subroutine LowRankRec is recursive. It takes as input the couple (A,r). Then: 

• if n < (m — r)^, it checks whether the algebraic set is empty or not; if "Dr = 0 it 
returns an empty list, otherwise it returns an error message; 

• if n = (m — r)^, it directly computes a rational parametrization of the projection of 
the hnite set Vr on (xi,..., Xn)] 

• if n > {m—rY, it computes a rational parametrization of the ideal (^(A o M, U, S, v)), 
where M, v are chosen randomly. Then it chooses a random value t G Q and recall 
recursively LowRankRec with input (A',r) where A' is obtained by setting Xi to t 
in A o M. 

The output is a rational parametrization whose real solutions contain a point for each 
connected component of 
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LowRankRec(74, r) 

Input: A n—variate linear matrix A{x) of size m, encoded by the m^(n+ 
1) rational entries of Aq, Ai,An, and an integer 0 < r < m — 1; 
Output: Either an error message if genericity assnmptions are not satis- 
hed, or a rational parametrization q = (go, <?i, • • •, Qn+i) ^ snch 

that for any connected component C C there exists to ^ ^(<?n+i) 

with (gi(to)/go(to), • • •, qn{to)/qo{to)) G C. 

Procedure: 

1. If n < (m — r)^ then retnrn an empty list; 

2. If n = (m — r)^ then retnrn Project(RatPar(/), (xi,..., Xn)); 

3. Choose randomly M G GL(n, Q), v G Q( 2 m-r)(m-r). 

4. P = lmage(Project(RatPar(/(A o M, U, S, v)), (xi,... ,x„)), M“^); 

5. Choose randomly t G Q, and dehne A = (Ao + tAi, A 2 ,..., An); 

6 . Q = Lift(LowRankRec(A,r),t); 

7. retnrn Union(Q, P). 


3.2 Correctness 

The correctness of the algorithm is proved in Theorem It is snpported by intermediate 
results that we discuss below. 

The hrst result states that Property G holds when A is generic in Also, it states 

that Property G is inherited by generic sections ft of the incidence variety. Its proof is 
contained in Section IH 

Proposition 2 Let S G GL(m — r, Q) and let U G 

1. There exists a non-empty Zariski-open set sZ C such that if A E D 

Qm2(n+i)^ (A, G, 5) satisfies Property G; 

2. If {A,U,S) satisfies Property G, there exists a non-empty Zariski open set PZ C C 
such that if t E AZ D Q, and A = (Aq + tAi) + X 2 A 2 + ■ ■ ■ + XnAn, then (A, U, S) 
satisfies G. 

The second result is about Lagrange system i{A o M, U, S,v). We prove that its set of 
solutions is hnite and that i satishes Property G 2 . The proof of this result is contained in 
Section [H 
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Proposition 3 Let S E GL(m — r, Q) and let U E be full-rank, and A such 

that {A, U, S) satisfies G. Let c = (2m — r){m — r). Then there exist non-empty Zariski 
open sets Y C and C GL(?7,, C) such that, ifvEYr\Q‘^ and M E fl the 

following holds: 

1. Z{A o M,Lf, S, v) is empty or finite and 1{A o M, U, S, v) satisfies Property G 2 ; 

2. the projection of Z{Ao M, U, S, v) on {x, y) contains the set of critical points of the 
restriction of tti to Vr ■ 

The following Proposition states that in generic coordinates the image of each connected 
component C of Vr A M" by each map Tifix) = {xi,... ,Xi) is a closed snbset of the real 
line M (for the Enclidean topology). In addition it states that the pre-images of valnes 
t G M lying on the boundary of 7ri(C) can be computed as projections of critical points of 
TTi restricted to Vr. 

Proposition 4 Let S E GL(m — r, Q). There exist non-empty Zariski open sets % C 
^ GL(n, C) such that if U E ^ n Q("*-px™- and M E ^2 G 
and if A is such that {A, U, S) satisfies G, and d = dAmVr, for any connected component 
C C Vr A M"", the following holds: 

1. for i = 1,... ,d, 7ii{M~^C) is closed; 

2. if t E M. lies on the boundary of 7ii{M~^C) then vrf^(t) fl M~^C is finite and there 
exists {x, y) E Vr{Ao M, U, S) such that 7ii{x, y) = t. 

Now, the previous results suggest that in order to prove correctness of the algorithm 
LowRank, both input data and parameters chosen internally must belong to pre-determined 
dense open subsets. We formalize this fact. 

If the input of LowRank is a (n + 1)— uple of matrices A = (Aq, Ai ,..., An), the computa¬ 
tion of the rational parametrization appearing in the procedure of LowRankRec is called 
exactly n — {m — r)^ -|- 1 times. We recall that we denoted by .s/ the non-empty Zariski 
open set computed in Proposition!^ Before the first call of LowRankRec, parameters U, S 
are chosen to define the algebraic set Vr{A,U, S) and no other choices of such matrices 
are done during the algorithm. 

Hypothesis Hi. Input A and random parameters U, S satisfy: 

• U is full-rank and S G GL(m — r, Q); 

• zi G 

One also has to ensure that U, S belong to the non-empty Zariski open sets defined in 
Proposition!^ and 1^ at each call of LowRankRec. Viceversa, parameter M strictly depends 
on the number of variables and so it is changed at each call. To summarize, the choices 
of random parameters done during LowRank can be stored in an array: 
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where the upperscript represents the number of variables. We also denote by 

^ 2 '^ and the non-empty Zariski open sets defined by Propositions and 4 
at the {n — j + 1)—th recursion call. 

Hypothesis H2. Random parameters (|^ satisfy: 


• U G r\"=im-rr ^ 

• S' G GL(m — r, Q); 

• G n n Q-^xj for j = (m — r)^,..., n; 

• G n Q(2m-r)(m-r) j _ _ ^^2^ . . . 

• G n Q for j = {m — r)^,..., n. 


Theorem 5 // Hi and H 2 hold, algorithm LowRank returns a rational parametrization 
whose set of solutions intersects each connected component ofVr flM". 

Proof : Suppose first that n < {m — r)^. Since Hi holds, then the variety Vr is empty. 
Hence fl M” = 0 and the algorithm returns a correct output. 

Thereafter, we proceed by induction on n > (m — r)^. 

If n = (m — r)^, since Hi holds. Proposition [^implies that Vr is finite. We deduce that the 
routine RatPar returns a rational parametrization of Vr and the routine Project returns a 
parametrization of "Dr- 

Let n > {m — r)^ and suppose that for any {n — 1)—variate linear matrix satisfying 
G, algorithm LowRank returns the expected output when Hi and H 2 hold. Let H be a 
n—variate linear matrix of size m, let r be an integer such that 0 < r < m — 1 and let 
C C ©rnM"" be a connected component. Let U, S be the random matrices chosen at Stepj^ 
of LowRank. Let M G GL(?7,, C) be the matrix chosen at Step|^of LowRankRec with input 
A and r, and consider the projection tti : (xi,..., Xn) —t xi restricted to Vr{A o M, U, S). 
Since Property Hi, H 2 hold, by Proposition]^ 7ri(M“^C) is a closed subset of the image, 
and so either 7ri(M“^C) = M or 7ri(M“^C) is a closed set with non-empty boundary. We 
claim that, in both cases, LowRank with input {A,r) returns a point which lies in the 
connected component M~^C. This is proved next. 

First case. Suppose hrst that 7 ri(M“^C) = M. In particular, for t G Q chosen at Step 
of LowRankRec with input H,r, the set 7 rj"^(t) intersects M~^C, so n M~^C 7 ^ 0. 

Let be the (n — 1)—variate linear matrix obtained from A o M hy substituting 

Xi = t. Remark that 7 rj"^(t) fl M~^C is the union of some connected components of 
the determinantal variety fl = {x E : rankHG-P < r}. Since Hi 

holds, U, S) satishes G ; we deduce that by the induction hypothesis the subroutine 

LowRankRec computes one point in each connected component of fl and so 

at least one point in M~^C. 

Second case. Suppose now that 7 ri(M“^C) 7 ^ M. By Proposition]^ 7 ri(M“^C) is closed. 
Since M~^C is connected, 7 ri(M“^C) is a closed interval, and since 'Ki{M~^C) 7 ^ M there 
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exists t in the boundary of 7ri{M~^C) such that 7ri(M“^C) C [t, +C)o) or 7ri{M~^C) C 
(—cxD,t]. Suppose without loss of generality that 7ii{M~^C) C [t,+cxD), so that t is the 
minimum value attained by tti on M~^C. 

Let X = {t,X 2 , ■ ■ ■ iXn) e M~^C. By Proposition |4| there exists y G such that 

{x,y) G Vr- We claim that there exists 2 ; G such that {x,y,z) lies in Z{A o 

M,U,S,v). Then, we conclude that the point x G M~^C appears among the solutions 
of the rational parametrization P obtained at Step of LowRankRec. Correction of the 
algorithm follows. Now we prove our claim. 

Let C' C Vr n pe the connected component such that M~^C' C M~^Vr H 

]^n+m(m-r) contains {x,y). We first prove that t = Tii{x,y) lies on the boundary of 
Indeed, suppose that there exists (x,y) G M~^C' such that 7ri(x, ^ < t. 
Since M~^C is connected, and since it is a connected component of a real algebraic variety, 
there exists a continuous semi-algebraic map r: [0,1] —)■ M~^C' with r(0) = {x^y) and 
r(l) = {x,y). Let be the map n^^x^y) = x. Then also vTa, o r is 

continuous and semi-algebraic (it is the composition of continuous semi-algebraic maps), 
and (tTx o r)(0) = x, {tIx o r)(l) = x. Since {tIx o 'r){0) G Vr for all 9 G [0,1], then 
X G M~^C. Since Tiiix) = 7ri(5;, ^ < f we obtain a contradiction. So iiiix^y) lies on the 
boundary of 

Since (A, U, S) satishes G, hence [A o M, U, S) satishes G, and by the Implicity Function 
Theorem one deduces that {x,y) is a critical point of tti restricted to Vr and that, by 
Proposition]^ there exists z G ^ Z{A o M,U, S,v), as 

claimed. □ 


3.3 Complexity analysis 


In this section we provide an analysis of the complexity of algorithm LowRank. We also 
give bounds for the number of complex solutions computed by LowRank. 

We suppose that f{A,U,S) satisfies Property G. We remark that this property can 
be checked, e.g. via Grobner bases techniques. Indeed, the dimension of sets Vr can 
be computed with Grobner bases. Further, using the Jacobian criterion [20l Theorem 
16.19], the regularity of Vr is easily checked by verifying that the complex solution set to 
f{A, U, S) and the maximal minors of its Jacobian matrix is empty. 


In order to bound the complexity of LowRank, it is essentially sufficient to bound the 
complexity of LowRankRec. This latter quantity mainly depends on the subroutine RatPar 


computing the rational parametrization, whose complexity is computed in Section |3.3.2 
We mainly rely on routines described in PI, which consists in a symbolic homotopy 
algorithm taking advantage of the sparsity structure of the input polynomial system. 


Finally, complexity bounds for the subroutines Project, Lift, Image and Union are provided 
in Section 3.3.3 and refer to results of 1551. 
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3.3.1 Bounds on the degree of the output of RatPar 


We consider the subroutine RatPar at the hrst recursion step of LowRank. Its input consists 
in either the generators f{A o M, U, S) of the incidence variety (if n = {m — r)^) or the 
Lagrange system i{A o M, U, S, v) (if n > (m — r)^). In both cases, we provide below in 
Proposition 1^ a bound on the degree of the rational parametrization returned by RatPar. 

We recall that if ..., are p groups of variables, and / G Q[x(^\ ..., we say 
that the multidegree of / is (di,..., dp) if its degree with respect to the group of variables 
is dj for j = 1,... ,p. 


Proposition 6 Let A be a n—variate m x m linear matrix, 0 < r < m — 1 and let U, S 
and M,v be respectively the parameters chosen at step of LowRank and at step [1 of 
LowRankRec. Suppose that Hi and H 2 hold. Then: 


1. if n = {m — r)^, the degree of the output of RatPar, with input f{A o M, U, S), is 
bounded from above by 6{m,n,r) = 6{m, (m — r)^,r) = 

2. if n > {m — r)^, the degree of the output of RatPar, with input 1{A o M, U, S), is 
bounded from above by 


5{m,n,r) = 


rr 


m{m — r)\ / n — 1 
n — k J \k + {m — r)^ — 1 


r(m — r) 
k 


with J^rn,n,r = {k '■ max{0 ,n — m{m — r)} <k < min{? 7 , — {m — r)‘^,r{m — n)}}. 


Proof of Assertion 1: If n = (m — r)^, since Hi holds, the dimension of the incidence 
variety Vr is zero. Consequently, the degree of the rational parametrization returned by 
RatPar is the degree of Vr- Since the entries of f{Ao M, U, S) are bilinear polynomials in 
X, y, one can compute the Multilinear Bezout bound (see |5S1 Chapter 11]) on this degree. 

From UY{y) — S one can eliminate {m—r)‘^ variables Pij (for example, those corresponding 
to the last m — r rows of Y{y)). Abusing notation, we denote by the same symbol / C 
Q[x, pi,i,..., yr,m-r\ the polynomial system obtained after this reduction. It is constituted 
by m{m—r) polynomials of multidegree bounded by (1,1) with respect to x = (xi,..., x„) 
and y (1/1,15 • • • 5 yr,m—r)- 


By |55l Proposition 11.1.1], degZ{f) is bounded by the sum of the coefficients of 
(s. + mod C s,]. 


Since n + rfm — r) = m{m — r), and since (s^, + is homogeneous of degree 

m{m — r), the aforementioned bound equals the coefficient of in the expansion 

of (s, + , that is . □ 


Proof of Assertion 2: In this case, the input of RatPar is the Lagrange system 
1{A o M,U, S,v). Let / be the reduced system defined in the proof of Assertion 1. 
We apply a similar reduction to 1{A o M,U, S,v). We introduce Lagrange multipliers 
z = [1, Z 25 • • • 5 Zm{m-r)] (we put Zi = 1 w.l.o.g., sffice i{A o M, U, S, v) is defined over the 
Zariski open set z ^ 0) and we consider polynomials {g, h) = z'Dif. Hence the new 
reduced system i = (/, g, h) is constituted by: 
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• m{m — r) polynomials of multidegree bounded by (1,1, 0); 

• n — 1 polynomials of multidegree bounded by (0,1,1); 

• r(m — r) polynomials of multidegree bounded by (1, 0,1). 


Moreover, by Proposition]^ Z(f,g,h) has dimension at most zero and {f,g,h) verifies 
G 2 . By [551 Proposition 11.1.1], degZ{f,g,h) is bounded by the sum of the coefficients 


of 


{s, + syr^^-^\sy + s,r-\s, + sj^^-^^ mod ^ Z[s,,Sy,S,]. 

As in the proof of Assertion 1, by homogeneity of the polynomial and by counting the 
degrees, the previous sum is given by the coefficient of the monomial 
in the expansion 

m(m—r) n—1 r(m—r) 


E E E 

i=0 j=0 k=0 


m[m — r ^ 
i 


n 


IN fT{jn 'f'}\ Yn{m—r)—i+j n—l—j+r{m—r)—k 


3 


k 




The coefficient is obtained by setting the equalities i + k = n, m{m — r)—i + j = r{m — r) 
and n — 1 — j + r{m — r) — k = m{m — r) — 1. These equalities imply i + k = n = 
j + k + {m — rY =j + k + i—j = i + k and consequently one deduces the claimed 
expression. □ 


Remark 7 If n = {m — rY and if the rank defect d = m — r is fixed, we have n = df is 
fixed and 5{m, n, r) = 5{m, df,m — d) = ^ 

Proposition 1^ implies straightforwardly the following estimate. 


Corollary 8 Suppose that the hypothesis of Proposition\^ are satisfied. Then LowRank 
returns a rational parametrization whose degree is less than or equal to 


/ m{m — r) 
\{m — rY 


mm{n,m^—r^} 

+ ^ <5(m,j,r). 

ji=(rn—r)2+l 


Proof : Since Hi holds, for n < {m — rY the algorithm returns the empty list. For m,j,r 
let J^m,j,r be the set of indices dehned in Proposition Observe that iFm,j,r = 0 if and 
only if j > m? — r^. Hence, the thesis is deduced straightforward from bounds given in 
Proposition □ 

One can also deduce the following bound on S{m,n,r). 

Lemma 9 For all m,n,r, with r < m — 1 , 6{m,n,r) < , 


Proof : This comes straightforwardly from the formula 



applied with a = n and b = m{m — r), and from the expression of d{m,n,r) computed in 
Proposition □ 
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3.3.2 Complexity of RatPar 


As announced in the preamble of the Section, our complexity model for RatPar is the 
symbolic homotopy algorithm in |^ . 

We suppose that n > [m — r)"^ and that the input of RatPar is the reduced Lagrange 
system i = 1{A o M,U,S,v) G Q[x,y, built in the proof of Assertion 2 of 

Proposition]^ One can easily build a second polynomial system i C Q[a;, y, z], such that: 

• the length of £ equals that of 

• for i = 1 ,..., n — 1 + w? — r^, the support of ii equals that of if, 

• the solutions of i are known. 

In fact, we remind that by construction, i contains three groups of quadratic polynomials 
in Q[a;,?/, 2 ;], of multidegree respectively bounded by (1,1,0), (0,1,1) and (1,0,1). We 
denote by Ai C Q[x, ?/], A 2 C Q[2/,2:] and A 3 C Q[a:,z] the supports of the three groups, 
so that for example Ai = { 1 , 0 :*, yj, Xiyj \1 <i <n,l <i < r{m—r)}, or, equivalently, Ai 
can be seen as the subset of made by the exponents of its monomials. Let £* be 

with support in Ai, 1 < i < m{m — r). Hence we generate two linear forms G Q[o:] and 
ki ,2 £ QM and we define ii{x,y) = ki^i{x)ki^ 2 {y)- We equivalently generate polynomials 
ii{y,z) = ki^i{y)ki^ 2 {z), m{m-r) + 1 < i < m{m-r)+n-l and ii{x,z) = ki^i{x)ki^ 2 {z), 
m{m — r)+n < i < n — l + m'^ — r^. We deduce straightforwardly that £ verifies the above 
properties, since its solutions can be computed by solving systems of linear equations. 

In llQ] , the authors build a homotopy path between i and i, such as 

ti + {1 — t)i G Q[x,y, z,t] (3) 

where t is a new variable. The system (j^ defines a 1—dimensional algebraic set, that 
is a curve. We deduce by [iQl Proposition 6.1] that, if the solutions of i are known, 
one can compute a rational parametrization of the solution set of system (j^ within 
0 ((n^AlogQ + h^^^)dd') arithmetic operations over Q, where: 

• n is the number of variables in i] 

• N = m{m — r)^Ai + (n — + r(m — is the cardinality); 

• Q = max*=i,2,3{||g|| : q e A*}; 

• d is the number of isolated solutions of i] 

• d' is the degree of the curve Z(ti + (1 — t)i)-, 

• ca is the exponent of matrix multiplication. 

Suppose the following preliminary lemma, whose proof is given below. 
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Lemma 10 Let Tm,,n,r be the set defined in Proposition^ and suppose J^m,n,r 7 ^ 0- Let 
6{m,n,r) be the bound defined in Proposition^ Then the degree of Z{ti + (1 — t)i) is in 

O ((n + — r^) min{?7,, m{m — r)} 6{m, n, r)) . 

Theorem 11 Letn > [m — rfi. Let A be an—variate mxm linear matrix, 0 < r < m —1 
and let M, U, S, v be the parameters chosen during the first recursive step of LowRank. Let 
5 = 5{m,n,r) be the bound defined in Proposition^ Then, RatPar returns a rational 
parametrization within 

0~ {fn + m^ — r'^y 5^) 

arithmetic operations. 


Proof : Following the notation introduced above, h = n — 1 + — r^. the bound for 

d is h and is given in Proposition and a bound for d' is given in Lemma [T^ and is in 
0~{n‘^6). Moreover, N G 0{nmr{m — ry), and hence N G 0{n^). The proof follows from 
HQl Proposition 6.1], since the maximum diameter of Ai, A 2 , A 3 is bounded above by n, 
that is Q < h. □ 

We exploit the multilinear structure of ti + {1 — t)i. By [55l 


Proof of Lemma 
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Proposition 11.1.1], degZ{ti + (1 — t)i) is bounded by the sum of the coefficients of 

(l{Sx^ Sy, Sz, St) = {Sx + Sy + + Sz + Si)"" ^(Sa; + Sz + 


modulo / = ^\sl) C Z[sx, Sy, Sz, St]. It is easy to check that q = 

qi + St{q2 + qs + qY + g with sf that divides g and 

qi = {Sx + Sy)^^^-^\sy + Sz)^-\sx + Sz)'-^^-^^ 

= m{m - r)st{sx + Sy)^^^-^^-ySy + Sz)^-\sx + 
q3 = {n- l)st{sx + Sy)^^^-^\sy + Sz)^-ySx + Sz)'’^^-^^ 
q, = r{m - r)st{sx + Sy)^^^-^\sy + Sz)^-\sx + 


and hence that q = qi + q 2 + qs + q 4 mod /. Below, we bound the contribution of 
qi,i = 1.. .4. The stated bound is given by the sum of the contributions and follows 
straightforwardly. 


Contributions of qi. The contribution of qi is the sum of its coefficients modulo the ideal 


n+l r{m-r)+l 


I' ={Sx . . 
with 6{m, n, r). 


Sz 


). This has been computed in Proposition 


and coincides 


The contribution of q 2 . Write q 2 = m{m — r)stq 2 with ^2 ^ '^[sx, Sy, Sz]. Consequently 
the contribution is given by the sum of the coefficients of q 2 , modulo /', multiplied by 
m{m — r). Now, observe that degg 2 = n — 2 + m^ —and that maxima powers admissible 
modulo P are s", Sy^™ \ Hence, three conhgurations give a contribution. 


(A) The coefficient of the monomial s” ^Sy*-™ ^ in ^ 2 , that is 

r(m—r) 


E 


A:=0 


m{m — r) — 1 


n — 1 


n — l — k J \k — 1 + (m — r) 


r{m — r) 
k 
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(B) The coefficient of the monomial ^ in q 2 , that is 

r(m—r) 

sb= 5: 


fc =0 


m{m — r) — 1 
n — k 


n — 1 

k — 1 + {m — r)^ 


r{m — r) 
k 


(C) The coefficient of the monomial s'^s^y'^ r)^m{m r) 2 

r(m—r) 

Ec= 


fc =0 


m(m — r) — 1 
n — k 


n — 1 

k — 2 + {m — r)2 


r{m — r) 
k 


So the contribution of q 2 equals m{m — 'r)(S^ + + Sc). 

One easily deduces that S^ < 5{m, n, r) and Sc < 5{m, n, r). Remember that we suppose 
J^m,n,r 7 ^ 0, that is (5(m, n, r) > 0. We claim that S^ < (1 + niin{n, m{m — r)}) 5{m, n, r). 
Consequently, we conclude that the contribution of q 2 is m{m — 'r)(S^ + S^ + Sc) € 
O {m{m — r) min{n, m{m — r)} 5(m, n, r)). 


Let us prove this claim. First, denote by 

Xi = max{0, n — m{m — r)} X 2 = min{r(m — r),n — {m — r)^} 

«! = max{0, n — 1 — m{m — r)} 0:2 = min{r(m — r), n — 1 — (m — r)^} 

the indices such that 5{m, n, r) sums over Xi < /c < X 2 and S^ over ai < k < a 2 . Remark 
that ai < xi and 0:2 < X 2 . Finally, denote by (f(k) the A;—th term in the sum dehning 
S^, and by j(k) the k—th term in the sum dehning S{m,n,r). 


For all indices k admissible for both S{m,n,r) and S^, that is for Xi < /c < 02 , one gets, 
by basic properties of binomial coefficients, that 

< '^(k)x(k) with Tf/c) = —^^ 

m{m — r)—n + l + k 

When k runs over all admissible indices, the rational function '^{k) is non-increasing 
monotone, and its maximum is attained in T(xi) and is bounded by mm{n,m{m — r)}. 
Three possible cases can hold: 



1. Xi = 0- Hence X 2 = n — {m — r)^, ai = 0 and a 2 = n — 1 — {m — r)^ = ai — 1. 
We deduce straightforwardly from the above discussion that < mm{n,m{m — 
r)}6{m, n, r); 

2. xi = n—m{m—r) and ai = n—l—m{m—r). We deduce that X 2 = ^2 = r{m—r) and 

that Sa = ‘d(^) ^ 9^(*^i) +niin{n, m{m — r)} 6{m, n, r) < (1 -|-min{n, m{m — 

r)}) (5(m, n, r); 

3. Xi = n — m{m — r) and ai = 0. Actually, either this case coincides with case 1 
(if n = m{m — r)) or we deduce that n = 1 + m{m — r), that is xi = X 2 = 1 and 
OL\= 0.2 = 0. By straightforward computations one deduces that Sa < S{m,n,r). 


The contribution of q^ and q^. Following exactly the same path as in the case of q 2 , one 
respectively deduces that the contribution of gs is in O {n mm{n,m{m — r)} 6{m,n,r)) 
and that of ^4 is in O {r{m — r) min{n, m{m — r)} S{m, n, r)). 

□ 
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3.3.3 Complexity of other minor subroutines 

For these complexity bounds, we refer to those given in [55l Lemma 10.1], [55l Lemma 
10.3], [55l Lemma 10.5] and [55l Lemma 10.6], from which they are obtained straightfor¬ 
wardly. 

Proposition 12 Let d(m,n,r) be the bound defined in Proposition^ At the first recur¬ 
sion step o/LowRankRec, the following holds: 

• the complexity of Project is in 0~ ((n -|- — r^fi n, r)Y); 

• the complexity of Lift is in 0~ {{n -L — r^) {S{m, n, r)fi); 

• the complexity of Image is in 0~ {{n -f m? — 5{m, n, r) + [n + — r'^fi); 

• the complexity of [)n\on is in 0~ {{n— r"^) {6{m,n,r)fi). 


4 Regularity of the incidence variety 


The goal of this section is to prove Proposition 
We start by proving Property Gi in Assertion 1 and 2. 

Proof of Property Gi in Assertion 1: We denote by dv C the algebraic set of 

mxm matrices of rank < r. By m Proposition 1.1] its singular locus is a^-i. For all A G 
C™' ("-+P, the set Vr is the intersection of of with the linear space C = Aq + (Ai, ..., Afi). 
By Bertini’s theorem (see [59l Theorem 17.16]), if C is generic, the following holds: 


sing Vr = sing(£ ft Or) = £ fl sing Or = CA fir-i = TAr-i- 


We conclude that there exists a non-empty Zariski open set C such that if 

A^ then A satisfies Property Gi. □ 

Proof of Property Gi in Assertion 2: We suppose that {A, U, S) satishes G. In 
particular A satisfies Property Gi. By Sard’s Lemma [55l Section 4.2], there exists a 
non-empty Zariski open set C C such that if t G then a point in Vr fl Z{x\ — t) 
is regular if and only if it is regular in Vr- Then, for t G .5^, the matrix obtained by 
instantiating xi to t in A satisfies Property Gi. □ 

We focus now on Property G 2 . We denote by ai^ij the entry of the matrix Ai at row i 
and column j, for £ = 0,1,..., n, z = 1, ..., m and j = 1,... , m. 

Our proof is essntially based on Thoms Algebraic Weak Transversality theorem [SSI Sec¬ 
tion 4.2]. We deduce, applying this result to given algebraic maps, the existence of a 
non-empty Zariski open set with the requested properties. 

Proof of Property G 2 in Assertion 1: Dehne 


/ ■ C”' X X 

{x,y,A) 


^m(m—r)+(m—r)^ 

f{A, U, S) 
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and, for a given A G define the induced map 

^ f(A,U,S). 


fA'. C" X 

{x,y) 


Suppose first that / ^(0) = 0. This is equivalent to saying that, for any A G 
Vr{A,U, S) = 0. By the Nullstellensatz [T71 Chap. 8], this implies that for any A G 
(j^m2(n+i)^ (/(A, f/, S)) = (1) which is a radical ideal, and so also (ft) = (1) for all f G C. 
In this case we conclude by defining ^2 = and ^ = C. 

Suppose now that /“^(O) 7 ^ 0. We prove that there exists a non-empty Zariski open set 
s:^2 C C™' C+ 1 ) such that if A G s^2-, the Jacobian matrix of f{A, U, S) has maximal rank 
at each point of /^^(O). This fact implies Property G 2 since we can apply the Jacobian 
Criterion [201 Theorem 16.19] to obtain that 

• the ideal generated by f{A, U, S) is radical; 

• the algebraic set defined by f{A, U, S) is either empty or smooth and equidimen- 
sional of co-dimension (2m — r){m — r). 

We claim that 0 is a regular value of /, i.e. at any point of the fiber /“^(O) the Jacobian 
matrix J associated to f{A, U, S) (with respect to variables x = (a:i,... ,Xn) and 
y = • • •) 2/m,m-r)) has maximal rank. By Thom’s Algebraic Weak Transversality 

theorem [HHl Section 4.2] we conclude that there exists a non-empty Zariski open set 
s^2. C such that, for every A G ^^ 2 , 0 is a regular value of the induced map /a, 

as claimed. We prove this claim in the sequel. 

Let (x, y, A) be in the fiber /“^(O) and J be the Jacobian matrix of /(A, U, S) evaluated 
at the point {x, y, A). We claim that there exists a maximal minor of J which is not zero 
at (x, y, A). In fact, we consider the submatrix of J obtained by: 

• the m{m — r) x m? block d{A{x)Y[y)) ■ j/dao^k,£ of derivatives of polynomial entries 
of A(^x)Y (y) with respect to the variables {ao,fcy | fc, £ = 1,..., m}; 

• the {m — X m{m — r) block d{UY(y) — S)/dy of derivatives of polynomial entries 
of UY{y) — S with respect to variables ?/; remark also that these polynomials do not 
depend on a. 

We remark that: the hrst block d(^A(^x)y)- j/dao^k,e, up to permuting rows and columns is 
a m{m — r) x block-diagonal matrix, with m blocks of size (m —r) x m on the diagonal 
all equal to Y(y)'; the second block d{UY(y) — s)/dy, up to re-ordering polynomials and 
variables, is a (m — r)^ x m{m — r) matrix with m — r blocks of size (m — r) x m on the 
diagonal all equal to U. 

Since {x,y,A) G /“^(O), U and Y(y) satisfy the matrix relation UY{y) = S and S is full 
rank. So U and Y(y) are full rank by the formula rank(f/y (y)) < min(rankf/, rank y (?/)). 
Moreover, polynomial entries of UY{y) — S' do not depend on parameters ao,ij- Hence 
we can extract a square non-singular submatrix of J of order m{m — r) + {m — r)^ = 
(2m — r)(m — r), proving that J has row-rank (2m — r){m — r). □ 
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Proof of Property G 2 in Assertion 2: Now, suppose that A E is given such 

that f(A, U, S) satisfies Property G 2 . Let 

TTi : Vr{A,U,S) -E C 
(x,y) ^ xi 

be the restriction to Vr of the projection on the first variable. Since / satisfies G 2 , 
Vr(A, M,s) is smooth and equidimensional of co-dimension (2m — r)(m — r). By Sard’s 
Lemma ( [551 Section 4.2]) the image of critical points of tti is included in a hypersurface 
of C. Let ^2 C C he the complement of this hypersurface. Then, if t G one of the 
following facts hold: 

• 7 r]"^(t) = 0. In this case = 0 and by the Nullstellensatz I(7r^^(t)) = (ft) = (1), 
which is a radical ideal; 

• Tilsit) 7 ^ 0 and for all {x,y) E {x,y) is not a critical point of the map 

TTi. So the Jacobian matrix of ft has full rank at each {x,y) E 2{ft), and by the 
Jacobian criterion ft defines a radical ideal and Z{ft) is smooth and equidimensional 
of co-dimension (2m — r)(m — r) -|- 1 . 


□ 

We conclude by proving that Vp is either empty or has dimension n — {m — p)^, and that 
this property is inherited after the recursive calls: this fact concludes Proposition 

Proof of dimension of Vp in Assertion 1: Let 0 < p < r. Let x denote the vector of 
m^ independent variables Xtj, I < i,j < m, and let X E be the matrix (x* j)i<i j<m- 

ByHa Proposition 1.1], the set 2^(minors(p J- 1, X)) C C™ is irreducible of co-dimension 
(m — pY dimension p(2m — p). Let be the entry of A^ at row i and column j. 
Let X = (xi,..., Xn) and let 

I = (minors(p 1,X)) -F {xtj - ao,ij - ai^tjXi -C Q[x,x]. 

The set 2^(minors(p J- 1,X)) C is irreducible of co-dimension [m — pY dimen¬ 

sion mf + n — {m — p)^ (in fact, variables Xi,..., do not appear). If linear forms in 

(Xjj - ao,ij - -are generic, then Z{I) C is empty or 

equidimensional of dimension n — {m — p)^ (Bertini’s theorem, see [5^1 Theorem 17.16]). 

Let tt: —?• C"" be the projection 7r(x, x) = x. Let V be an irreducible component of 

Z{I). Then V has dimension n — {m—pY. For x G vr(l/), the fiber 7r“^(x) is finite. By the 
Theorem on the Dimension of Fibers [56l Sect. 6.3, Theorem 7], dim 7r(l/) = dim V = 
n — {m — pY. We conclude that there exists a non-empty Zariski open set C C™ 
such that if A G then Vp has dimension n — {m — p)^. We conclude by defining 

Proof of dimension of Vp in Assertion 2: Now, suppose that A E is 

such that Vp has co-dimension (m — p)^ for some 0 < p < r. Hence, by Bertini’s 
theorem ([^ Theorem 17.16]) there exists C C such that if t G D Q, and 

A = (Aq J- tAi) + X 2 A 2 J- ■ • • -t- XnAn, hence Vp = {x E | ranky4(x) < p} has 

co-dimension {m — pY or is empty. We conclude by defining AA^ = 
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The proof of Proposition is now immediate by defining fl - 0 ^ fl and = 

n ^2 n ^3- 


5 Dimension of Lagrange systems 

The goal of this Section is to prove Proposition!^ This will be done in Section [5^ Before 
that, we give a local description of the incidence variety and of the Lagrange system and 
we provide intermediate results. 


5.1 Local description of the incidence variety 


Let A = Aq + xiAi + ... + XnAn be a n—variate m x m linear matrix with coefficients in 
Q, and let r < m — 1. Here we give a local description of the incidence variety Vr{A, U, S). 
From now on, for / G Q[x], we denote by Q[a:]/ the local ring of Q[a:] at /. We recall 
that the polynomial system defining Vr is given by /(H, t/. S'), which contains the entries 
of A{x)Y{y) and UY(y) — S. For p < r, let N be the upper-left p x p submatrix of A, so 


that 


A 


N Q] 
P' R ) 


(4) 


with Q G P' G and R G p). The next Lemma 

computes the equations of Vr in the local ring Q[x, j/]det Af- 


Lemma 13 Let A, N, Q, P, R be as above, and U, S be any full-rank matrices. Then 
there exist {qi,j}i<i<p,i<j<ra-r C Q[a:]detAf and {qlj}i<i,j<m-p C Q[a:]detw such that the 
constructible set Vr H {{x,y) : det A^(a:) 7^ 0} is defined by the equations 

yi,j - - • • • - qi,m-pym,j = 0 l<i<p, 

+ • • • + (l'i,m-pymj = 0 l<i<m-p,l<j<m-r 
UY{y) -S = R 


Proof : We denote by Y^^'^ and Y^"^^ the submatrices of Y(y) containing respectively the 
first p rows and the last m—p rows. We also use the block-division of H as in Q. We claim 
that in Q[x,|/]det 7 v the m(m — r) equations A{x)Y(y) = 0 can be read as the m{m — r) 
equations: 


/ IpYW + iV-iQy(2) X ^ 

V E(iV)F( 2 ) J - u 


where S(iV) = R — P’N~^Q is the Schur complement of N in A. Renaming the entries of 
N~^Q and Ti{N) concludes the proof. To prove the claim, remark that since det iV 7 ^ 0, 
A{x)Y(y) = 0 if and only if 


Ip 0 \ ( N-^ 0 

-P' T Jin T 


N Q 
P' R 


Y{y) = 0. 


□ 


23 



5.2 The rank at a critical point 

Given A, N, P, Q, R, S(A^) as above, let 


A = 


4 iV-iQ \ 

0 E(N) ) • 


Lemma 13 implies that the equations of Vr in the open set {{x,y) : det iV 4 0} can be 
rewritten as AY(y) = 0 and UY(y) — S = 0: the polynomial entries of the above ex¬ 
pressions are elements of the local ring Q[a:]det n- Now, from the first group of relations 
A[x)Y{y) = 0 one eliminates variables {yi,j}i<i<p,i<j<m-r-i which can be expressed as 
polynomial functions of x and {yi^j}p+\<i<m,i<j<m-r- That is, using the notations intro¬ 
duced in Lemma 1^, we can express the entries of Y^^'^ as polynomials in x and in the 
entries of Y^‘^\ 


Now, consider relations UY(y) — S = 0 where the entries of Y^^'^ have been eliminated. 
This is a linear system in the entries of with coefficients in Q[a;]det 7 v- Since S is 
full-rank, then U is full-rank and hence UY{y) — S = 0 consists of (m — r)^ independent 
relations. Finally one can eliminate (m — r)^ among the (m — p){m — r) entries of 
(suppose the first (m — r) rows) and re-write E{N)Y^‘^'^ = 0 as (m — p){m — r) relations 
in X and in the last (r — p){m — r) entries of Y^^\ 

Call F this polynomial system, consider Lagrange multipliers z = (zi,..., Z(^rn-p)(m-r)) 
and the polynomial system 


{gi, ...,gn)= z'D^F - {wi,.. .,Wn). 

The solutions to the above polynomial system contain the critical points of the projection 
■Kw restricted to Vr H {(x, y) : det N 4 0}. The next Lemma shows that, when w is generic 
in C”, the solutions to the Lagrange systems project on points of Fr with the expected 
maximal rank r. 

Lemma 14 LetA,U,S be as above and suppose that (A, r) satisfies G 3 . Letp <r — l and 
let g = {gi ,..., gn) be the polynomial system defined above. Then there exists a non-empty 
Zariski open set W <Z such that if w E W then g = 0 has no solution. 

Proof : Let C C be the constructible set defined by = 0 and 

by det iV 4 0 and rank A{x) = p, and let C be its Zariski closure. Let : (x, p, z,w) ^ x 
be the projection on the first n variables. The image 'XxiC) is included in Fp C F^ 
and hence, since {A,r) satisfies G 3 , it has dimension at most n — [m — p)^. Moreover 
this projection is dense in Fp. The fiber of tIx over a generic point x E Fp is the graph 
of the polynomial function w = z'DxF, and so it has co-dimension n and dimension 
{r — p){m — r) + {m — p){m — r) = {m — r){m + r — 2p). By the Theorem of the Dimension 
of Fibers [5^ Sect. 6.3, Theorem 7] one deduces that the dimension of C (and hence of 
C) is at most n — (m — p)^ -|- (m — r) (m -|- r — 2 p) = n — {r —pfi. Since p < r — 1 then C 
has dimension at most n — 1. We deduce that the projection of C onto the space C"" of 
tc is a constructible set of dimension at most n — 1 , and it is included in a hypersurface 
H E FA. Defining W = F'^\H ends the proof. □ 
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5.3 An intermediate lemma 


We consider the incidence variety Vr = Vr(A, S) and the restriction of the projection 
Tiw : (x, y) —)■ w'x to Vr, with w G C"'. Under the hypothesis that (A, [/, S) satishes G, the 
set Vr is either empty or smooth and eqnidimensional of co-dimension c := (2m—r)(m—r). 
Hence the set of critical points of the restriction of to Vr is the projection on the (x, y)- 
space of the solntion set of the system 

where 2; = {zi,..., Zc, Zc+i) 7^ 0 . Now, for A G at any solution {x,y,z) of this 
polynomial system, Zc+i 7^ 0 . Moreover, with the hypothesis w ^ 0, also (zi, ..., Zc) 0- 
So one introduces a linear combination ^2^=1 ~ 1 with v G Q'^, and we consider the 

system 

C 

/ = 0 , g = 0, h = 0, '^ViZi-l = Q. ( 5 ) 

i=l 



The polynomial system ([^ has n + c- 1 -m(m — r) -|-1 polynomials and n-l-c-t-m(m — r)-|-l 
variables. We denote by yVw{A,U, S,v) C +1 its zero set. We prove in next 

Lemma that when w is generic, the set yVw{A,U, S,v) is finite and regular, and that it 
encodes the critical points of restricted to Vr. The proof exploits the local description 


of Vr given in Section 5.1 


Lemma 15 Let (H, U, S) satisfy G. There exist non-empty Zariski open sets Y (Z and 
W <Z LIT such that if v and w G W, the following holds: 

1. the set yVw{A, U, S,v) is finite and the Jacobian matrix of ([^ has maximal rank at 
each point of Ww{A, U, S, v); 

2. the projection of Ww{A,U, S,v) on the {x,y)-space contains the critical points of 
TTw ■ {x, y) —t w'x restricted to Vr(A, U, S). 


Proof of Assertion 1 of Lemma 15: We first show that the Lagrange system (|^ 


can be re-written in a local form when we consider the local description of the incidence 
variety Vr as above. Let W C C" be the set defined by Lemma 14l and w G W. Then one 


has that all solutions (x, y, z) to the system ([^ are such that rank A(x) = r. Hence, there 
exists ar xr submatrix N of A{x) such that det A 7 ^ 0. We prove below that there exist 
non-empty Zariski open sets yfj C C'^ and C C” such that for v G Yff and w G #jv, 
the statement of Assertion ( 1 ) holds locally. Hence, to retrieve the global property, it is 
sufficient to define y (resp. W) as the hnite intersection of sets yfj (resp. Wn), where 
N varies in the collection of r x r submatrices of A. 


We use the block-division of matrix A as in (|^ with p = r and without loss of generality 
one can assume to work in the open set det A 7 ^ 0, with A the upper-left r x r sub matrix 
of A. We deduce by Lemma that the local equations of Vr are 

= -N-^QY^‘^\ S(A)V(2) = 0, f/(by(i) + f/(2)y(2) _ 
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where is the row-subdivision of the matrix Y(y) as in Lemma 13 and is 


the correspondent column-subdivision of U. From the hrst and third groups of equations 
one obtains that S = U^^\-N-^QY^‘^^) + f/(2)y(2) = (_f/(i)iv-iQ + f/(2))y(2). Since S 
is full-rank, then ^nd —U^^'^N~^Q + are non-singular, and so: 

• the second group of equations can be re-written as S(A^) = 0; 

• the third group of equations can be re-written as Y^"^^ = 

The entries of S(A^) in the local ring Q[a:]detv are exactly the (m — r)^ minors of A{x) 
obtained as determinants of the (r 4-1) x (r -1- 1) submatrices of A{x) containing N (see 
for example the proof of [55l Proposition 3.2.7]). Since (A, U, S) satisfies G, A satisfies Gi, 
and the Jacobian Dx[T,{N)]ij of the vector of entries of S(iV) has full-rank at each point 
X such that rank74(a;) = r. 

We call f = (/(,..., /') the local equations represented by the entries of Y^^'> + N~^QY^'^\ 
S(A^) and The Jacobian matrix of /' has the form 


Df = {DJ' Dyf) = 


We consider the polynomials 


0(m—r)2 xm(m—r) 


^r{m—r) 

0 




idlJ • • • ) 9ni ^1) • • • ) ^m(m—r)) ('^1; • • • ) ^c+l) ^ 


DJ' Dyf 
Wi... Wn 0 


Polynomials in h' = (h),..., h'^(m-r)) relations Zi = 0, for i = (m — rf + 1,..., c, 

and can be eliminated together with variables Zi,i = (m — rf -|- 1,... ,c. So the local 
equations of the Lagrange system ([^ are: 

(m—r)^ 

f = 0, g' = 0, ^ ViZi -1 = 0, 

i=l 

for a given v G This is a square system with n -|- c + 1 polynomials and n -|- c -f 1 

variables. Now, consider the map 


p : C’"+‘=+^ X C" X 

{x,y,z,w,v) 


and its section map 


{x,y,z) 


(J^n+c+l 

(f'.g'.T.tYv.Zi-i) 


(j^n+c+1 


for given v G and w G C"". If p ^(0) = 0, then for all v,w, = 0, and the 


claim is proved by taking F' = T)(r = and W = Wn = W (see Lemma 14). 
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Suppose now p“^(0) 7 ^ 0 and let {x,y, z,w,v) G p“^(0). We claim that the Jacobian 
matrix of p at {x,y, z,w,v) has maximal rank. Hence, 0 is a regular value for p and by 
Thom’s Weak Transversality Theorem [SSI Section 4.2] there exist C and 

Wn C C” non-empty Zariski open sets, such that if n G and w G then 0 is a 
regular value for This implies that, by the Jacobian criterion, the set W^(H, U, S, v) fl 
{{x,y,z) I detiV(a;) 7 ^ 0} is empty or zero-dimensional, and that the claim is true with 
X (Cc-(m-r) ^ -vy-picP ig algo non-empty and Zariski open. We prove below this 
claim by exhibiting a non-singular submatrix of Dp. 

We remark that, since (H, U, S) satishes G, the Jacobian matrix Df has maximal rank at 
{x,y). Moreover, Zc+i 7 ^ 0 and by the relation '^ViZi — 1 = 0 there exists 1 < £ < (m —r)^ 
such that zg 7 ^ 0. Then we consider the submatrix of Dp obtained by isolating: 


• the non-singular submatrix of Df'] 

• the derivatives of gi,... ,gn with respect to wi,... ,Wn, giving the identity block 

Idjii 

• the derivative of '^ViZi — 1 with respect to vg. 


The previous blocks isolate a matrix of size (n -|- c -f 1) x (n -f c -|- 1) whose determinant 
does not vanish at {x, y, x, w,v). □ 


Proof of Assertion 2 of Lemma 15 


This proof is similar to that of |3Zl Lemma 15]. 


Suppose hrst that yVw{A,U, S,v) = 0 for all tc G C"" and v G C'’. Fix w G C”, 
{x,y) G Vr{A,U, S) and suppose that {x,y) is a critical point of Tiy^ restricted to Vr- 
Then, since Vr is equidimensional, there exists x 7 ^ 0 such that {x, y, z) verihes the equa¬ 
tions z'Df = [tc,0]. Since ^ 7 ^ 0, there exists n G C'^ such that v'z = 1. So we conclude 
that {x,y, [ 2 :,!]) G Ww{A,U, S,v), which is a contradiction. Hence crit(7r^„, Vr) = 0 and 
Assertion 2 is proved. 


Suppose now that {A,U,S) satisfy G and that Z{p) is non-empty. Suppose that w E W 
(the set dehned in the proof of Assertion 1). By [55l Sect. 3.2], crit(7rr„, Vr) is the image 
of the projection 7ix,y on x,y of the constructible set: 


S = {{x,y,z) : f = g = h = 0, z^O} 

where f^g, h have been dehned in (|^. One can easily prove, by means of Thom’s Weak 
Transversality theorem, that S has dimension at most 1. Moreover, for each {x,y) in 
the hber 7ify{x,y) has dimension 1, since if {x,y,z) G 7if^{x,y), then {x,y,Xz) G 
7T~y{x,y) for all A 7 ^ 0. By the Theorem on the Dimension of Fibers [56l Sect. 6.3, 
Theorem 7], we deduce that 7ix,y{S) is hnite. Fix now {x,y) G 7ix^y{S) and let Y(x ,y) 
be the non-empty Zariski open set such that if n G Y(x,y), the hyperplane YhviZi — 1 = 0 
intersects transversely 7r“y(a;,j/). Let Y' C C'’ be the set dehned in the proof of Assertion 
1. By dehning 

y = r fl 

(x,y)eTT^^y{S) 

one concludes the proof. Indeed, ^ is a hnite intersection of non-empty Zariski open sets. 

□ 
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5.4 Proof of Proposition 


We finally use Lemma to show that, up to a generic change of variables, the set of 
critical points is hnite. 

Proof of Proposition 1^ Let C GL(n, C) be the set of M G GL(n, C) such that the 
first row w' of M~^ lies in the set W given in Lemma this set is non-empty and Zariski 
open since the entries of M~^ are rational functions of the entries of M. Let Y <Z 
be the non-empty Zariski open set given by Lemma 15 and let v ^ "V . We denote by 
e; = (l,0,...,0)eQ". 

Remark that for any M G the following identity holds: 


( Df{AoM,U,S)\^( Df{A,U,S)oM\( M 0 \ 

V e'i0---0 ) w'0---0 /V^ Imim-r) ) 


where Df{A, U, S) o M means that at all entries of Df{A, U, S) we substitute x ^ M x. 
We conclude that the set of solutions of the system 


f{A, U, S) = 0, (^ 1 ,..., z,)Df{A, U, S) + zc+i{w', 0) = 0, v'z-l = 0 (6) 


is the image by the map 

/ X \ / M~^ 0 0 \ / X \ 

1 / ^ 0 Im{m -r) 0 2 / • 

\ z J \ 0 0 Ic+i J \ z J 

of the set S of solutions of the system 

f{AoM,U,S) = 0, {z,,...,z,)Df{AoM,U,S) + z,+,{e[,0) = 0, uG - 1 = 0. (7) 

Now, let 71 be the projection that forgets the last coordinate Zc+i- Remark that vr(5) = 
Z{A o M,U,S,v) and that tt is a bijection. Moreover, it is an isomorphism of affine 
algebraic varieties, since if {x, y, z) G S, then its 2 ;c-i-i-coordinate is obtained by evaluating 
a polynomial at (x, y,zi... Zc). 

Thus, Assertion 1 of Lemma [T^ implies that: 


• S and vr(5) = Z{A o M,U, S, v) are hnite; 

• the Jacobian matrix associated to ([^ has maximal rank at any point of S. 


Since vr(5) = Z{Ao M, U, S, v) and that vr is an isomorphism. Assertion 1 of Proposition 
m follows. 


Assertion (2) is a straightforward consequence of Assertion 2 of Lemma 15 


□ 


6 Closure properties 

The goal of this section is to prove Proposition We use notation of 123 Section 5], 
which we recall below. 
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Notations Let ^ C be an algebraic variety of dimension d. The z— eqnidimensional 
component of Z is denoted by i = 0,... ,d. We denote by y'iZ) the nnion of the 

following sets: 


• Qq{Z) U ■ ■ ■ U Qd-i{Z) 

• the set sing(f2(i(Z)) of singnlar points of 

and by ^(vrj, Z) the Zariski closnre of the nnion of the following sets: 


• Qo{2) U ■ ■ ■ U Qi-i{Z); 

• the nnion for r > z of the sets crit(7rj, reg(f2r(^))) of critical points of the restriction 
of TTj to the regular locus of Qr{Z). 

For M G GL(zz,C) we recursively dehne the collection of algebraic sets 
as follows: 

• Od{M-^Z) = M-^Z; 

. 0,{M-^Z) = U'^(7r,+i,0,+i(M-iZ)) U^(7r,+i,M-iZ) for z = 

0,...,d- 1. 

We recall the two following properties dehned in [371 Section 5]. 


Property P{Z). Let ^ C C"" be an algebraic set of dimension d. We say that M E 
GL„(C) satishes P(Z) when for all z = 0,1,..., d 

1. Oi{M~^Z) has dimension < z; 

2. Oi{M~^Z) is in Noether position with respect to xi,..., x*. 


Property Q{Z). Let Z be an algebraic set of dimension d and 1 < z < d. We say that 
Qj(^) holds if for any connected component C of the boundary of 7rj(C) is contained 

in TTi{Oi_i{Z) n C). We say that Q(^) holds if Qi{Z) holds for all 1 < z < d. 

In [27] the authors proved that given any algebraic variety Z of dimension d, Property 
P{Z) holds generically in GL(rz, C) (Proposition 17) and that if M G GL(rz,C) satishes 
P{Z), then Q{M~^Z) holds (Proposition 18). We use these results in the following proof 
of Proposition 

Proof of Assertion 1 of Proposition Let ^2 C GL(n, C) be the non-empty 
Zariski open set computed by m Proposition 17] for Z = T>r- One obtains that any 
M G ^2 verihes P(T’r)- Remark that since M G GL(zz, C) there is a natural bijective 
correspondence between the set of connected components of Vr O M” and the ones of 
OM” given by C GG M~^C. Fix a connected component M~^C C M~^Vr flM"' and 
consider the projection tt* restricted to M~^Vr flM"'. Since M G ^ 2 , by m Proposition 
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18] the boundary of 'Ki{M ^C) is contained in 7ii{Oi-i{M ^Vr) flM ^C) and in particular 
in 7ii{M~^C). This implies that 7rj(M“^C) is closed and so Assertion 1. □ 

Proof of Assertion 2 of Proposition Let us prove Assertion 2. Let M E ^2 and 
let t e M lie in the boundary of ni{M~^C). By [371 Lemma 19] the set fl M~^C is 

hnite. Fix x E 7rh^(t) H M~^C, and let p < r be the rank of A{x). Since x is hxed, the 
polynomial system y 1 —)■ f{A, U, S) parametrized by U, S is linear in y and can be written 
in the form 

(] ]) Y(y) = ]) (8) 

Let B be the matrix on the right side of (|^. This system is equivalent to m — r linear 
systems of equations whose unknowns are the columns of Y[y) and whose constant terms 
are the columns of B. Now, each of these systems has a solution if and only if, by RouchA 
Capelli Theorem [13], each column of B lies in the space spanned by the columns of A{x), 
that is if and only if the following equality holds: 


rank 


A{x) 

U 


rank 


A(x) r) \ 

U S )' 


Denote by fhe number of rows of U that do not lie in the space spanned by the 
rows of A{x) and by fhe number of rows of \U \ S'] that do not he in the space 

spanned by the rows of [A{x) \ 0]. Since S is full-rank, one necessarily deduces that for 
all U, rx,u < 'fx,u,s = m — r (in fact, no line of S is the zero vector). Moreover, since x 
is hxed, then there exists a non-empty Zariski open set ^c,x such that HUE ^c,x then 
= m — r and 


rank 


A{x) 

U 


^ =P + rx,u 


p + m — r 


and 


rank 


u s 


^ =p + rx,u,s 


p + m — r, 


so that the system has at least one solution. One concludes by dehning 


^ = n n 

CcX’rnR’^ 3::G7r^^(£)nC 


which is non-empty and Zariski open by the finiteness of the number of connected com- 
ponents of Vr 0 M"" and of the set 7rj|“^(t) fl C. □ 

7 Practical Experiments 

This section reports on practical experiments made with a hrst implementation of our 
algorithm. Note that for computing rational parametrizations, we use Grobner bases and 
change of ordering algorithms [251 EH- Our experiments are done using the C library 
FGb, developed by J.-G. Faugere [23] and interfaced with Maple. The implementation 
will be freely released as a Maple library and made available at 
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http://homepages.laas.fr/henrion/software 


We start by comparing our implementation with implementations of general algorithms 
based on the critical point method in RAGlib [52]. Next, we comment the behaviour 
of our algorithm on special examples that are well-known by the research community 
working on linear matrices. 


7.1 Comparison with RAGlib 


We have generated randomly linear matrices for various values of m and n and run our 
implementation for different values of r. By randomness of rational numbers we mean 
that we generate couples of integers chosen with uniform distribution in a hxed interval. 
Clearly, this would imply that the set of inputs is hnite (hence, it is not a Zariski dense 
set). On the other hand, this does not affect genericity and also the correctness of the 
algorithm since the requested properties can be checked before its execution. 

Our implementation is written in Maple. As said above, we use the Grobner engine 
FGb for computing in practice rational parametrizations. All computations have been 
done on an Intel(R) Xeon(R) CPU A7540@2.00GHz 256 Gb of RAM. We report in Table 
numerical data of our tests. For any choice ofm, 2<r<m — 1 and n, we generate a 
random dense linear matrix A and we let LowRank run with input {A,r). 


(m, r, n) 

PPC 

LowRank 

deg 

maxdeg 

(m, r, n) 

PPC 

LowRank 

deg 

maxdeg 

(3,2,2) 

0.2 

6 

9 

6 

(5,2,3) 

0.9 

0.5 

0 

0 

(3,2,3) 

0.3 

7.5 

21 

12 

(5,2,4) 

1 

0.5 

0 

0 

(3,2,4) 

0.9 

9.5 

33 

12 

(5,2,5) 

1.6 

0.5 

0 

0 

(3,2,5) 

5.1 

13.5 

39 

12 

(5,2,6) 

3 

0.6 

0 

0 

(3,2,6) 

15.5 

15 

39 

12 

(5,2,7) 

4.2 

0.7 

0 

0 

(3,2,7) 

31 

16.5 

39 

12 

(5,2,8) 

8 

0.7 

0 

0 

(3,2,8) 

109 

18 

39 

12 

(5,2,9) 

OO 

903 

175 

175 

(3,2,9) 

230 

20 

39 

12 

(5,3,2) 

0.4 

0.5 

0 

0 

(4,2,2) 

0.2 

0.5 

0 

0 

(5,3,3) 

0.5 

0.5 

0 

0 

(4,2,3) 

0.3 

0.5 

0 

0 

(5,3,4) 

43 

22 

50 

50 

(4,2,4) 

2.2 

2.5 

20 

20 

(5,3,5) 

OD 

5963 

350 

300 

(4,2,5) 

12.2 

26 

100 

80 

(5,4,2) 

0.5 

125 

25 

20 

(4,2,6) 

OO 

593 

276 

176 

(5,4,3) 

10 

167 

105 

80 

(4,2,7) 

OO 

6684 

532 

256 

(5,4,4) 

OO 

561 

325 

220 

(4,2,8) 

OO 

42868 

818 

286 

(5,4,5) 

oo 

5574 

755 

430 

(4,2,9) 

OO 

120801 

1074 

286 

(6,3,3) 

4 

1 

0 

0 

(4,3,3) 

1 

8 

52 

36 

(6,3,4) 

140 

1 

0 

0 

(4,3,4) 

590 

18 

120 

68 

(6,3,5) 

OO 

1 

0 

0 

(4,3,5) 

(OO 

56 

204 

84 

(6,3,6) 

OO 

2 

0 

0 

(4,3,6) 

OO 

114 

264 

84 

(6,3,7) 

oo 

2 

0 

0 

(4,3,7) 

oo 

124 

284 

84 

(6,3,8) 

oo 

2 

0 

0 

(4,3,8) 

oo 

124 

284 

84 

(6,4,2) 

0.6 

40 

0 

0 

(4,3,9) 

oo 

295 

284 

84 

(6,4,3) 

1 

64 

0 

0 

(4,3,10) 

oo 

303 

284 

84 

(6,4,4) 

341 

300 

105 

105 

(4,3,11) 

oo 

377 

284 

84 

(6,5,3) 

95 

276 

186 

150 

(5,2,2) 

0.6 

0.5 

0 

0 

(6,5,4) 

(OO 

8643 

726 

540 


Table 1: Timings and degrees for dense linear matrices 


We compare our timings (reported in column “LowRank”) with the function PointsPer- 
Components (column “PPC”) of the real algebraic geometry library RAGlib, implemented 
by the third author [52] • The symbol oo means that the no result has been returned after 
4 days of computation. 


31 









We make the following remarks about Table [T] 


1. We first observe that our algorithm is most of the time faster than RAG lib and it 
allows to tackle examples that are out of reach of RAGlib. 


2. The growth in terms of timings with respect to n seems to respect the correspondent 
growth in terms of degrees of output parametrizations ; in particular note that we 
have established that for r and m hxed, the sum of the degrees of parametrizations 
we need to compute stabilizes when n grows. This is observed in practice of course 
and is reflected in our timings compared to those of RAGlib. 


3. Accordingly to the related Multilinear Bezout Bounds computed in section 3.3.1 , the 
degrees of rational parametrizations stabilize when n grows, since when n > m? — 
and the input is generic, LowRank does not compute critical points at hrst calls. 
This fact is remarkable, since: 


• a natural geometric invariant associated to its degree as complex alge¬ 
braic set, does not depend on the dimension n of the affine section (one can 
prove easily that generically this degree is given by Thom-Porteous-Giambelli 
formula, cf. [21 Ch. II,§4]); 

• an algebraic invariant naturally associated to the output-size (the degree) is 
constant in n, coherently with the abovementioned geometric invariant. 


Finally, we give a hnal remark on potential a posteriori verihcation of the correctness of the 
output of LowRank. Deciding whether a hnite set, encoded by a rational parametrization, 
meets every connected component of a given real algebraic set, is a hard problem, far 
from being solved, both from a theoretical and computational viewpoint. As far as the 
authors know, there are no symbolic or numerical algorithms able to perform this task. 
Also, producing such a certihcate seems to be hard to imagine, but this was not among 
the goals of this paper. Anyway LowRank is able to produce a certihed correct output 
when regularity assumptions (that can be checked a priori) holds. 


7.2 Examples 


In this last section, we consider some examples of linear matrices coming from the liter¬ 
ature, and we test the behavior of LowRank. We consider examples of symmetric linear 
matrices since, as observed in Section [h2 the main motivation for solving the real root 
hnding problem is to obtain dedicated algorithms for spectrahedra and semidehnite pro¬ 
gramming. 


Example 16 (The Cayley cubic) We consider the 3x3 linear matrix 


A{x) = 
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The real trace of the complex determinantal variety V 2 = {x G \ rank A{x) < 2'^ is 
shown in Figure^ The convex region {x G | A{x) F 0} is the Cayley spectrahedron. 
We run LowRank with input {A,r) with r = 2 and r = 1 (the case r = 0 is trivial since 

A[x) is always non-zero and hence Vq is empty). In both cases, the algorithm first verifies 

that the genericity assumptions are satisfied. 

Let us first analyze the case r = 2. LowRank runs 3 recursive steps. Its output is a 
rational parametrization of degree 14 with 12 real solutions and 2 complex solutions. We 
give below details of each recursive call of LowRankRec. At the first, at a rational 

parametrization of degree 5 is returned, with the following 5 real solutions: 

r/l\ /1\ /-1\ f ~^\ ( 18.285118452 \ ] 

< 1,-1, 1 , -1 , 164.322822823 ) . 

[\1/ V-1/ V-1/ Vw \ 4.552268485 / J 

The coordinates of the fifth point are approximated to 9 certified digits and such approxi¬ 
mation can be computed by isolating the coordinates in intervals of rational numbers as: 

^ 1- 21081306277346124211 21081306277346754459 ] ^ i q 1 Q/l 

"''I ^ 1. 1152921504606846976 ’ 1152921504606846976 J ~ iO.ZOOllO^OZ 
^ r 5920353629066611305 23681414516266799197 ] ^ 

■^2 L 36028797018963968 ’ 144115188075855872 J ~ lU^.OZZOZZOZO 

^ [ 10496816461511385723 2624204115377866059 ] ^ a 
"''S t [ 2305843009213693952 ’ 576460752303423488 J ~ 

Remark that it also computes the 4 singular points ofV 2 , where the rank of A is 1. At the 
second (resp. third) recursive call, it returns a rational parametrization of degree 6 (resp. 
of degree 3) with 4 (resp. 3) real solutions. 

In the case r = 1, step [2] of LowRankRec returns a rational parametrization of degree 4 
which encodes the 4 singular points 0 /P 2 H that is Vi fl At the second and third 
recursions, LowRankRec return empty lists. 

We finally remark that the above results are typical, in the sense that the 4 singular points 
contained in Vi fl are always computed at the first recursion step, both in case r = 2 
and r = 1. Conversely, the coordinates of the other real solutions depend on the choice 
of random parameters (while their number is constant). Moreover, all computations end 
after few seconds (< 5 sec.). 


Example 17 Let 


/ Oi 

Xi 

X 2 

X3 \ 

Xi 

02 

X3 

X4 

X2 

3^3 

O3 

X5 

\ Xs 

X4 

X5 

O4 / 


where x = {xi,X 2 ,x^,Xi,xfi) are variables and ( 01 , 02 , 03 , 04 ) are parameters. We first 
instantiate (01,02,03,04) to (1,1,1,1) and in this case we obtain that with input {A,2) 
and (^,3), the genericity assumptions reguested by LowRank are not satisfied. For r = 1, 
these assumptions are satisfied and the algorithm returns a degree 4 parametrization, with 
4 real solutions encoding Pi fl 


Subseguently, we let (oi, 02, 03, 04) vary randomly in For all random instantiations, we 
observe that the inputs {A, 3), [A, 2) and {A, 1) verify the genericity assumptions, and that 
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the degrees of the rational parametrizations returned at each recursion step are constant, 
while the number of real solutions changes with parameters. We summarize our results in 

TablelB 



r = 3 

r = 2 

r = 1 

partial degrees 

[12 24 24 12 4] 

[12 20 8 0 0 ] 

[0 0 0 0 0 ] 

total degree 

76 

40 

0 

time (s) 

768 

21.5 

4.6 


Table 2: Degrees and timings for Example 17 with generic parameters 


Example 18 (The pillow) Let 

xi 0 xi \ 

1 0:2 0 

X 2 I X 3 

0 X 3 1 J 

The spectrahedron S' = {x G | A{x) ^ 0} is known as the pillow, see also 
Section 5 . 1 . 1 ], It is pictured in Figure^ with the help of the software povray (http: 
//WWW. povray. org ) implementing the raytracing algorithm. 


A{x) = 


/ 1 

Xi 

0 

V Xi 



Fignre 2: The pillow and its algebraic bonndary 

The Zariski closure of its boundary is the real trace of the complex hypersurface defined 
by the vanishing of 

det A{x) = 1 — xl — xl — 2x1 + ^1^3 ~ 2 xlx 2 X 3 + x\x 2 . 

As clear from Figure\^ the determinantal hypersurface consists in four branches arising 
from the convex set S. 

The boundary of S contains 4 singular points of the determinantal hypersurface, where 
A{x) has rank 2. Their coordinates can be found by computing a Grobner basis of the ideal 
generated by the 3x3 minors of A, which is {2x\ — 1, 20:3 — 1,X2 + X 3 }. In particular, 
these four points are contained in the hyperplane X 2 = —X 3 . 
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We tested LowRank with input {A, 2). We obtain that at the first recursion, at step a 
rational parametrization q = (q'o, ? 2 , ? 3 , ^ 4 ) of degree 4 (with only real roots) is com¬ 
puted. By isolating the 4 real roots of as in Example\T^ one gets the following rational 
approximations of the singular points: 


Xi e 
X2 e 
X3 e 


r 6521908912666475339 13043817825332644843 ] ^ 
1. 9223372036854775808’ 18446744073709551616 J ~ 


-V2/2 


[ 26087635650665343561 6521908912666428733 
1-36893488147419103232 ’ 9223372036854775808 


] ^ V2/2 


I 6521908912666412349 13043817825332731855 ] ^ 
1. 9223372036854775808’ 18446744073709551616-1 ~ 


-V2/2. 


As for Example \T^ we observe a typical output in terms of the degree of the rational 
parametrizations and the number of real solutions. Details are given in Table\^ 



r = 3 

r = 2 

r = 1 

partial degrees 

[6 8 4] 

[4 0 0] 

[0 0 0] 

total degree 

18 

4 

0 

real solutions 

14 

4 

0 

time (s) 

< 5 

< 5 

< 5 


Table 3: Degrees and timings for Example 18 for the pillow 
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